X<-read.table("data/prostata.tab",header=TRUE)
X<-X[-1,]
p<-ncol(X)
Y<-as.numeric(X[,1])
Y[Y==2]<-1
Y[Y==3]<-0
X<-matrix(as.double(as.matrix(X[,-1])),nrow=length(Y), ncol=p-1)
v<-apply(X,2,sd)
m<-apply(X,2,mean)
X<-t(apply(t(apply(X,1,"-",m)),1,"/",v))
n<-nrow(X)
X<-cbind(1,X)

